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ABSTRACT 

We have performed a search for old globular clusters (GC) using archival F606W and F814IU 
HST/ACS images of 19 Magellanic-type dwarf Irregular (dlrr) galaxies. Those dlrrs reside in nearby 
(2 - 8Mpc) associations of only dwarf galaxies. All dlrrs have absolute magnitudes fainter than or 
equal to the SMC (My = —16.2 mag). We detect in total 50 GC candidates in 13 dlrrs, of which 
37 have (V — I) colors consistent with "blue" (old, metal-poor) GCs (bGC). The luminosity function 
(LF) of the bGC candidates in our sample shows a turnover magnitude of My = —7.41 ± 0.22 mag, 
consistent with other galaxy types. The width of the LF is a = 1.79 ± 0.31 which is typical for dlrrs, 
but broader than the typical width in massive galaxies. The half-light radii and ellipticities of the 
GCs in our sample (r\ ~ 3.3 pc , e ~ 0.1) are similar to those of old GCs in the Magellanic Clouds and 
to those of "Old Halo" (OH) GCs in our Galaxy, but not as extended and spherical as the Galactic 
"Young Halo" (YH) GCs (r~h ~ 7.7 pc , e ~ 0.06). The e distribution shows a turnover rather than a 
power law as observed for the Galactic GCs. This might suggest that GCs in dlrrs are kinematically 
young and not fully relaxed yet. The present-day specific frequencies of GCs (Sn) in the galaxies in 
our sample span a broad range: 0.3 < Sn < 11. Assuming a dissipationless age fading of the galaxy 
light, the Sn values would increase by a factor of ~ 2.5 to 16, comparable with values for early-type 
dwarfs (dE/dSphs). A bright central GC candidate, similar to nuclear clusters of dEs, is observed in 
one of our dlrrs: NGC 1959. This nuclear GC has luminosity, color, and structural parameters similar 
to that of uj Cen and M 54, suggesting that the latter might have their origin in the central regions of 
similar Galactic building blocks as the dlrrs in this study. A comparison between properties of bGCs 
and Galactic YH GCs, suspected to have originated from similar dlrrs, is performed. 
Subject headings: globular clusters: general — galaxies: star clusters — galaxies: dwarf — galaxies: 
irregular 



1. INTRODUCTION 

Hierarchical structure formation models envision the 
assembly of present-day massive galaxies via numer- 
ous minor mergers/accretions of smaller galactic entities. 
Present day dwarf galaxies have masses that are similar 
to those predicted for the proto-galactic fragments (e.g. 
Purcell et al. 2007), which were later incorporated into 
more massive galaxies (e.g. Searle & Zinn 1978). The 
prediction of the hierarchical growth of major galaxies 
through merging of many dwarf-sized fragments at early 
times is manifested by the steepening of the faint-end 
slope of the galaxy luminosity function with redshift (e.g. 
Ryan et al. 2007; Khochfar et al. 2007). Furthermore, 
dwarf irregular galaxies are found to be one of the most 
abundant galaxy types in the high-z Universe (e.g. El- 
lis 1997; Stiavelli et al. 2004), and thus are regarded as 
building blocks of massive galaxies. As nearby dlrrs are 
less evolved systems than nearby massive galaxies, they 
might represent the most "pristine" dwarfs that arc likely 
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similar to the fragments incorporated into massive halos 
and therefore hold important clues as to how galaxies 
formed. 

Old globular clusters (GC) are among the first objects 
to form in the early Universe. As luminous agglomera- 
tions of coeval stars with homogeneous metal abundances 
they represent unique tools for tracing the main galax- 
ies' star formation episodes. Therefore, their properties, 
as distinct entities and as a system, reflect the physical 
conditions at the time of their early formation. This will 
provide us with constraints on the galaxies' early assem- 
bly, since dlrrs may well contribute to the assembly of the 
rich GC systems (GCSs) of the most massive galaxies. 

Numerous observations of the most massive elliptical 
and spiral galaxies in various environments have shown 
that they contain exceedingly rich populations of old GCs 
(e.g. Kundu et al. 1999; Larsen et al. 2001; Dirsch et al. 
2003; Harris et al. 2006; Tamura et al. 2006). Extensive 
studies were performed in the last decade to understand 
how such populous GCSs were assembled (e.g. Ashman & 
Zepf 1992; Zepf & Ashman 1993; Hilker et al. 1999; Goud- 
frooij et al. 2003; Chandar et al. 2004; Puzia et al. 2004; 
Rhode et al. 2005; Goudfrooij et al. 2007). These studies 
led to two important discoveries: i) the bimodal metal- 
licity/color distribution of the GCs (e.g. Ashman & Zepf 
1992; Neilsen & Tsvetanov 1999; Gebhardt & Kissler- 
Patig 1999; Puzia et al. 1999; Kundu & Whitmore 2001) 
and ii) the presence of young/intermediate age massive 
star clusters in merging, starburst and irregular galaxies 
(e.g. Whitmore & Schweizer 1995; Goudfrooij et al. 2001; 
Puzia et al. 2002; Goudfrooij et al. 2004). As a result, 
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three major galaxy/GCS assembly scenarios have been 
proposed to explain these findings. The first scenario 
is the hierarchical build-up of massive galaxies through 
merging and accretion of pre-galactic dwarf-sized gas 
fragments (Searle & Zinn 1978) in which the metal-poor 
GCs form in situ while the metal-rich GCs originate from 
a second major star formation event (Forbes et al. 1997) 
from infalling gaseous fragments, e.g. the mini-mergers 
at high redshifts (Beasley et al. 2002). The spiral-spiral 
dissipativc galaxy merger scenario (Schweizer 1987; Ash- 
man & Zepf 1992) assumes that the metal-poor GCs were 
formed early in the progenitor galaxies, while metal-rich 
GCs formed during the major merger events. Finally, 
the accretion scenario incorporates the classical mono- 
lithic collapse picture in which the galaxies and their 
GCs form in the initial starburst (Pipino et al. 2007)and 
the later accretion of smaller dwarf-sized systems con- 
tributed to the assembly of massive galaxies GCSs (Cote 
et al. 1998; Hilker et al. 1999; Cote et al. 2002; Beasley 
et al. 2002; Kravtsov & Gnedin 2005). All three scenarios 
are not mutually exclusive and actually bear a number 
of similarities; they vary mostly in the amount of gas 
involved. For extensive discussion on the strengths and 
weaknesses of the various scenarios we refer the reader 
to the reviews by Ashman & Zepf (1998); Kissler-Patig 
(2000); van den Bergh (2000); West et al. (2004); Brodie 
& Strader (2006). 

Dwarf irregular galaxies play a key role in each 
GCS/galaxy formation scenario. As dlrrs host mainly 
old metal-poor, blue GCs (bGCs), the GCSs of massive 
galaxies, and mainly their bGCs, might be a mixture of 
clusters formed in situ and GCs accreted with the dwarfs 
through hierarchical merging. Thus, it is very important 
to study dlrr galaxies to obtain a proper understanding 
of their GCS properties and what might be their contri- 
bution to the blue GC population of massive galaxies, 
which has so far not yet been possible due to the lack of 
large enough samples. 

Analyzing GCs' metallicities and horizontal branch 
(HB) morphologies, Zinn (1993) proposed a sub- 
classification of the Milky Way (MW) GCS into three 
populations - BD (bulge/disk), OH (old halo) and YH 
(young halo), where BD are all metal- rich ([Fc/H] > 
—0.8) clusters, and the metal-poor, OH and YH, are di- 
vided by their horizontal branch index, i.e. AHB > —0.3 
and AHB< —0.3, respectively. With the idea that the 
YHs might be of external origin, first Zinn (1993) and 
later Mackey & Gilmore (2004) and Mackey & van den 
Bergh (2005) showed that in many respects the prop- 
erties (colors, luminosities, spatial distributions, struc- 
tural parameters, HB morphologies etc.) of the MW YH 
GCs are similar to the properties of GCs in dwarf galax- 
ies (LMC, SMC, Fornax, Sagittarius). Recently, Lee 
et al. (2007) showed that GCs with extended horizontal 
branches (EHBs) are decoupled in their orbital dynam- 
ics and mass from the rest of the GCs and confirmed the 
Zinn (1993) conclusion that the YH clusters are domi- 
nated by random motions and share the hot kinematic 
properties of EHB Galactic GC population supporting 
their probable external origin. This, and evidences of 
tidal streams around the MW and M31 (e.g. Ibata et al. 
2001; Martin et al. 2004; Kalirai et al. 2006) further sup- 
ported the view by Searle & Zinn (1978) that a fraction 
of the GCs were formed in isolated low-luminous satellite 
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building blocks, which were later accreted. 

The present study focuses on the GC properties (lumi- 
nosity, color, structural parameters, spatial distribution 
and specific frequencies) of dlrr galaxies in groups and as- 
sociations of only dwarf galaxies. The observational data, 
its reduction, GC selection and photometry are described 
in Sect. 2. The results are presented and summarised in 
Sect. 3, where the properties of GC candidates are dis- 
cussed and compared with the Galactic YH GCs. Ad- 
ditionally, wc address the question whether these dlrrs 
can be the progenitors of the present day dE/dSph by 
comparing their GC specific frequencies (SW) with those 
of early-type dwarfs in clusters (e.g. Miller et al. 1998; 
Seth et al. 2004; Forbes 2005; Strader et al. 2006; Miller 
& Lotz 2007). 

2. OBSERVATIONS 

2.1. HST/ACS data 

The current study uses HST archival data of 20 
Magellanic-type dwarf irregular (dlrr) galaxies located 
in nearby (d = 3 — 8Mpc) loose groups (associations) of 
dwarf galaxies (Tully et al. 2006). They were observed 
in the F606W and F814W filters with HST/ACS dur- 
ing Cycle 13, program GO-10210 (PI: B. Tully). The 
split into two non-dithered exposures was designed to 
reach ~ 2 mag beyond the Tip of the Red Giant Branch 
(TRGB) luminosity at Mi = -4.05, {M v = -2.5) (Da 
Costa & Armandroff 1990; Lee et al. 1993) in order to de- 
rive the TRGB distances to these groups/associations of 
(dlrr) galaxies. Accurate distances to these galaxies are 
published in Tully et al. (2006), which provides us with 
an excellent tool to study the properties of the GCSs of 
these dwarfs. 

2.2. Initial photometry 

The retrieved archival images were processed with the 
standard HST/ACS archival pipeline. In order to im- 
prove the object detection and initial photometry, the 
diffuse galaxy light was modeled and subtracted by con- 
volving the images with ring aperture median kernel of 
41 pixel radius. Then, on the residual images we used 
the IRAF 6 DAOPHOT/DAOFIND routine to detect 
objects with > 4-er above the background. Photom- 
etry for all detected sources was performed with the 
DAOPHOT/PHOT routine in apertures of 2,3,5 and 
10 pixel radius. To convert the instrumental magnitudes 
to the ST magnitude system an aperture correction from 
10 pixel aperture to infinity was applied, using the zero 
points and aperture correction values from the ACS pho- 
tometric calibration by Sirianni et al. (2005). The magni- 
tudes derived this way are labelled in the following with 
the subscript 'lOtot'. These magnitudes were only used 
for the GC selection, while the final GC magnitudes were 
determined as described in Section 2.5. 

2.3. Aperture corrections 

The typical size of a GC, containing half of its light, 
is about 3pc. Therefore, the GCs in the galaxies in our 
sample are expected to be resolved on ACS images, and 

6 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — Constant aperture correction is color and magnitude 
dependent. In panel a) are shown the magnitude difference de- 
rived from the individual curves of growth for each GC candidate 
(F606Wtot magnitudes) versus the magnitudes simply corrected 
with a constant value from 10 pixels to infinity (F606Wiotot )• 
Analogously, in panel b) we show the same relation for the ob- 
jects' colors. The thick line is the least-square linear fit to the data 
and 7 is the slope of the fit. It is seen that, even when exclud- 
ing the faintest objects, which are also the reddest (F606Wtot - 
F814Wtot) > 0.1), a simple aperture correction is magnitude and 
color dependent. 



typically their light extends beyond 10 pixels, which can 
introduce systematic errors in the aperture correction. 
We explore how this error affects our GC selection cri- 
teria. In Figure 1 are shown the differences between the 
corrected (F606Wiotot; see above) aperture magnitudes 
(Fig. 1 a) and colors (Fig. 1 b) and the "true" magnitudes 
(F606Wtot) derived from curves of growth for each GC 
(see Sect. 2.5). The Sirianni et al. (2005) aperture cor- 
rections are computed for point sources, however, in the 
case of resolved GCs this correction will be insufficient 
and introduce an offset in the GCs' magnitudes. Figure 1 
shows that this effect appears to be magnitude and color 
dependent. The dashed lines indicate the zero offset, the 
thick lines a linear fit to the data and 7 is the slope of the 
fit. The observed magnitude dependence can be under- 
stood in the sense that the brighter the object the more 
extended it is. Thus, the correction is underestimated 
for brighter GCs while it is overestimated for the fainter 
ones, when a constant aperture correction was applied 
to all GC candidates. The color dependence of the aper- 
ture correction is even more significant (Fig. lb). Even 
when the reddest ((F606Wi tot-F814Wi tot) > 0.1) ob- 
jects are excluded as likely background galaxies, the color 
dependence does not change significantly. Hence, a con- 
stant aperture correction could affect the clusters' colors 
by up to 0.1 mag. This result has significant implications 
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Fig. 2. — The initial GC selection was based on SSP evolution- 
ary models in the ACS/WFC filter system for different metallicities 
(Anders & Fritze-v. Alvensleben 2003). Namely, objects with col- 
ors -0.4 <(F606W-F814W)< 0.2, i.e. objects older than > 4 Gyr 
for the whole range of metallicities. With rectangles are shown the 
typical color ranges for Galactic metal-poor (MP) and metal-rich 
(MR) GCs. 

for studies dealing with similar quality data. 

2.4. GC Candidate Selection 

Our goal is to study the old GC population in these 
dlrrs. Thus, our GC candidates selection criteria tar- 
gets GCs with ages older than ~ 4 Gyr for a range of 
metallicities as provided by the GALEV Simple Stellar 
Population (SSP) models in the ACS filter system (An- 
ders & Fritze-v. Alvensleben 2003). Therefore, we re- 
strict the color selection limits to —0.4 < (F606W — 
F814W) S tmag < 0.15. As seen in Fig. 2, the limit at 
(F606W-F814W)stmag ~ -0.4 mag (corresponding to 
Johnson/Cousins V-lw 0.7) is 0.1 STMAG bluer than 
the most metal-poor GALEV model. With rectangles in 
Fig. 2 are shown the typical color ranges for metal-poor 
(MP) and metal-rich (MR) Milky Way GCs. Compar- 
ison with Bruzual & Chariot (2003) SSP models show 
that evolutionary tracks with 0.02 x Zq and 0.005 x Zq 
have V — I colors literally indistinguishable for ages 
older than 4 Gyr. Therefore, our color limit at (F606W- 
F814W)stmag > —0.4 selects GC candidates with ages 
older than 4 Gyr and (V — I) > 0.7 mag. Considering the 
weak dependence of (V — I) colors on metallicity, our GC 
candidates might be as metal-poor as 0.005 x Zq. 

Since our GC selection is based on effects mentioned 
in Sect. 2.3 as well. Given the small dispersion in the 
faint magnitude bins seen in Fig. 1 a, the effect of the 
aperture correction on the initial GC selection is negli- 
gible. Therefore, we adopt a faint magnitude selection 
limit at the TRGB, i.e. M v = -2.5 (Mi = -4.05) 
(Da Costa & Armandroff 1990; Lee et al. 1993), which is 
~ 5 mag fainter than the typical GC luminosity function 
turnover magnitude at My, to = —7.4 mag (e.g. Harris 
1996, 2001). This absolute magnitude was converted to 
apparent magnitudes using the distance modulus for each 
galaxy as derived by Tully et al. (2006). 

To correct for foreground Galactic extinction we have 
used the E(B — V) values towards each galaxy from the 
Schlegel et al. (1998) dust maps 7 . We used the Galactic 
extinction laws by Cardclli et al. (1989) (their eqs. 2 and 
3) to compute the values for the ACS filters. 

7 http:/ /nedwww. ipac.caltech.edu/ 
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Fig. 3. — Summary of the initial GC selection process for UGC 685. In Panel a) are shown the color and magnitude selections. They 
are indicated by a rectangular region (dashed box) for all objects obeying the shape selection criteria (imexam e < 0.15; FWHMj?go6w — 
FWHMj?8i4w). In panels b) and c) are shown the objects selected this way which were analyzed against their FWHM and light concen- 
tration (both m2 — m3 and ni2 — ms). In panel d) are shown the two concentration parameters against each other. The objects within the 
dashed box were selected as GC candidates. Panel e): the selected GC candidates (circled) are shown with different symbol types. The 
candidates which after visual inspection were regarded as GCs are labeled. See Sect. 2.4 for detailed explanation. 
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As mentioned above, due to the deep imaging and 
high ACS spatial resolution, the GCs in our images show 
profiles significantly more extended than point sources 
(stars). In order to include this knowledge in the GC 
selection we used IMEXAM to obtain initial information 
on the profiles of all detected objects as first approxima- 
tion for the GC selection. We considered as GC can- 
didates round objects (FWHMp 60 ew ~FWHM f8 i Wl 
ellipticity< 0.15) with 2 < FWHM < 9 pix (the measured 
stellar FWHM was ~ 1.7 to 1.9 pix). The imexam e and 
FWHMs were measured with fixed r = 5 pix aperture ra- 
dius and Moffat index j3 = 2.5, typical for stellar profiles. 
Hence, their absolute values cannot be taken as final re- 
sult but rather as upper limits for the initial selection 
criteria. The final cllipticitics and half-light radii (e, r^) 
were measured with ISHAPE (Sect. 2.6). As seen in Fig. 15 
and 17, the derived e values for final GC candidates reach 
0.3, a typical upper value observed for old GCs in the 
Magellanic Clouds and our Galaxy (Fig. 16, 17). There- 
fore, the adopted for initial selection IMEXAM e < 0.15 
within r = 5 pix, should not have introduced an artificial 
bias in our selection, as can be expected for such resolved 
objects. An additional check, with an increased e selec- 
tion limit, proved to include large number of contamina- 
tion due to blends and nuclei of background galaxies. 

The selection process is illustrated in Figure 3, with 
UGC 685 as an example. In Fig. 3 a) we show the color- 
magnitude diagram (CMD) for all UGC 685 objects sat- 
isfying the aforementioned size selection criteria. The 
dashed box indicates the color and magnitude selection 
limits (see above). The observed dispersion in our CMD 
is simply due to crowding effects leading to the uncer- 
tainties in the 10 pix photometry. Hence, stars affected 
by their neighbors in crowded regions will show 'wrong' 
colors and increased color dispersion. 

As the GCs are expected to show concentrations dif- 
ferent to stars and background galaxies, we incorporate 
this information by calculating the difference between the 
2, 3 and 5 pixel aperture photometry. In terms of con- 
tamination, the ni2— m.3 index is the more robust star- 
GC-background galaxy discriminator than the ni2— nis 
index. This is due to the increase of contaminating 
light from resolved Galactic stars with increasing aper- 
ture size. The GCs occupy a specific locus in the FWHM 
vs. m2-ni3 or ni2— 1115 plane, typical for resolved objects. 
In Fig. 3b) and c) we show the concentration 1112- m.3, 
1112— 1x15 indices versus the measured FWHM^si4vk for 
all color, magnitude and shape selected objects (imexam 
e < 0.15, FWHM F6CW -FWHM f81w ). The dividing 
line around m2-ni3 ~ 0.4 and ni2— 1TL5 __ 0.75 mag sep- 
arates stars from GC candidates. However, the bulges 
of resolved big galaxies and barely resolved compact 
(FWHM< 9 pix) background galaxies fall in the same re- 
gion as the GCs (see also Puzia et al. 2004). In Fig. 3d) 
the two concentration indices are plotted against each 
other. The dotted box shows the region occupied by the 
GC candidates. It can be seen that if we were to select 
GC candidates from the ni2— ni5 index alone, we would 
have introduced more contamination (shown with blue 
dots) from blended sources within the galaxy. The com- 
pleteness and contamination issue is discussed in more 
detail in Sect. 2.7. 

In the case of UGC 685, the automated GC selection 
procedure returned 13 objects, which were visually in- 



spected. 5 obvious background galaxies and 2 blends 
were removed leaving 6 GC candidates which are shown 
with numbered IDs in Fig. 3e) (see also Fig. 6). Four 
of the five objects with 1112-3 — 0.4 and rri2_5 ~ 0.7 
in Fig. 3 e) were visually inspected and turned out to be 
blended sources in the galaxy body, one was a clear back- 
ground galaxy. The same was the case in the rest of the 
studied galaxies, therefore we adopted the lower limit for 
concentration selection at the quoted values above. The 
same technique was applied to all our sample galaxies. 
The final GC candidate sample comprises 60 objects in 
13 out of the 19 examined dlrr galaxies. 

After the presented selection process, the contamina- 
tion in our GC sample by unresolved background galaxies 
cannot be ruled out entirely Completeness, however, is 
not an issue due to the depth of the images, with 90% 
point source limiting magnitude at My = —2.5 mag. 
Therefore, we arc confident that we do not introduce 
any unaccounted artificial bias in our GC selection. A 
similar selection of GC candidates in nearby LSB dwarf 
galaxies was recently confirmed with 96% success rate by 
radial velocities (Puzia & Sharina 2008). 

2.5. Final Globular Cluster Candidate Photometry 

Accurate photometry of the total light of resolved GCs, 
including virtually all cluster light, was based on individ- 
ual curves of growth which account for the different GC 
sizes due to the varying distances to the host galaxy. 
However, in images where a significant part of the host 
galaxy stellar component is resolved a special treatment 
of the contaminating sources (galaxy stars, complex stel- 
lar regions, foreground Galactic stars, background galax- 
ies), found in the photometric apertures, is needed. In 
order to eliminate the contaminating sources and build 
curves of growth up to 50 pixel aperture radius, we used 
the SExtractor (Bertin & Arnouts 1996) package. Mea- 
suring the GCs' magnitudes with SExtractor turned out 
to be an unstable task leading to very uncertain measure- 
ments especially for objects in crowded regions. Instead, 
we used the SExtractor option "CHECKIMAGE.TYPE 
-OBJECTS" to detect and subtract iteratively in two 
steps contaminating objects having different profiles. For 
their detection two different convolution filters were used, 
provided within SExtractor, i.e. gauss_5.0_9x9.conv and 
gauss_l.5_3x3.conv. 

The iterative procedure is illustrated in Figure 4 where 
two cases with strong contamination are shown. In the 
lower row the brightest star on the right was masked 
out. The two iterative steps were designed to detect and 
subtract only sources with the given filter profiles from 
the images and replace them with the background value 
as evaluated by a 60-pixcl background mesh. This al- 
lowed us to adjust the detection thresholds such that (as 
a first step) the GC with the associated pixels within 
0.5cr of the background (thus practically loosing no light 
from the GC) was subtracted from the original image. 
The image of this step, containing only contaminating 
objects, (Fig. 4 b) was then subtracted from the origi- 
nal (Fig. 4 a), thus removing only the contaminants. The 
second iteration used the narrower convolution filter and 
the object detection was setup such that the GC was 
detected and subtracted from the result of the first it- 
eration (Fig. 4 c). Then this image was subtracted from 
the result of the first iteration, thus removing the last 
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Fig. 4. — An illustration of the reduction process of iterative cleaning of the photometric apertures around the GC candidates in two 
extreme cases with strong contamination by resolved galactic stars and a nearby bright foreground star. As discussed in the text, two 
cleaning iterations (columns b) and c) ) using SExtractor were involved to remove the underlying galaxy light and contaminating sources 
within the photometric apertures. The CCD image noise characteristics were restored with MKNOISE (column d). Aperture photometry 
and curves of growth were build from the images reduced this way. 



contaminants (see Figure 4 d). Since SExtractor replaces 
the detected sources with a constant background value, 
we used the ACS CCD characteristics from the image 
header and restored the detector noise in the image with 
the IRAF MKNOISE routine. This step was mainly per- 
formed to restore the proper error budget in the pho- 
tometry. The circles in Fig. 4 mark some of the aperture 
radii used in our curve of growth photometry. 

In order to build individual curves of growth, an aper- 
ture photometry from 5 to 50 pixels with a step size of 
5 pixels was performed with the DAOPHOT/PHOT 
task on the final images (free of contaminants). Fifth- 
order polynomial fits to the curves of growth were used 
to derive the magnitudes for each GC. In order to com- 
pare our GC magnitudes with those of other GC studies 
we iteratively converted the instrumental STMAG mag- 
nitudes from ACS/WFC filters to the Johnson/Cousins 
V and / magnitudes using the synthetic transformation 
coefficients computed by Sirianni et al. (2005). 

2.6. Globular Cluster Sizes 

With the spatial resolution (0.05 arcsec/pixel) of the 
ACS/WFC the GCs in these nearby (2 to 8 Mpc) dlrrs 
are clearly resolved. This enables one to measure their 
structural parameters with relatively high precision. 

The GC half-light radii (r^) were estimated using 
the ISHAPE algorithm (Larsen 1999), which models the 
source as an analytic function convolved with a (model) 
point-spread function (PSF). For the PSF model we used 
the TinyTim software package 8 . It properly takes into 
account the field-dependent WFC aberration, filter pass- 
band effects, charge diffusion variations and varying pixel 
area due to the significant field distortion in the ACS field 
of view (Krist & Hook 2004) . We have used the GCs co- 
ordinates to create ten times sub-sampled model PSFs at 

8 http:/ /www. stsci.edu/software/tinytim/tinytim. html 



the position of each GC. When sub-sampling is enabled, 
TinyTim doesn't convolve the PSF with the charge dif- 
fusion kernel (CDK), but provides its values in the PSF 
image header. As it is position and wavelength depen- 
dent we generated an individual CDK for each GC, which 
was later used with ishape in the profile fitting process. 

Assuming a King (1962) model for the intrinsic radial 
luminosity profile of the clusters, the clusters ^60611^ 
images were modeled with all concentration parameters 
C = r t /r c = 5, 15, 30, 100 provided by ISHAPE. An 8 pixel 
fitting radius was used. The FWHM of the best y 2 model 
was adopted for the objects. To test the stability of the 
measured FWHMs we varied the fitting radius in the 
range 4 to lOpix corresponding from <~ 6 to 15 pc. No 
significant trend and practically identical FWHMs were 
returned by ishape. In order to convert from FWHM 
to rh we used the conversion factors tabulated for KING 
models with different C by Larsen (2006). However, this 
relation is valid for circularly symmetrical profiles while 
the output of ishape is the FWHM along the major axis 
(w y ). In order not to overestimate the rh by simply as- 
suming its measured rh,w value, a correction consider- 
ing the ratio of the FWHMs along the minor and major 
axes (w x /wy) should be applied. Therefore we adopted 
for the "true" effective radius the geometric mean of the 
FWHMs along the two axes, i.e. 



Hi = 'V 



c/Wy 



(1) 



where w x /w y is the ishape output of the minor/major 
axes FWHMs ratio. The comparison between values 
derived this way and the numerically derived correction 
(eq. 1 1 in Larsen 2006) resulted in absolutely identical 
values. The analysis and interpretation of the GCs sizes 
derived as described above is presented in Sect. 3.3. 

2.7. Completeness and Foreground/Background 
Contamination 
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Fig. 5. — Apparent V-band luminosity (panel a) and color (panel 
b) distribution of the expected Hubble UDF background contam- 
ination (filled histograms). The dashed-dotted line, in panel a) 
shows the luminosity distribution of all GC candidates in our sam- 
ple while with thick (blue) and dashed (red) histograms show the 
candidates divided by their color into blue and red GCs, respec- 
tively (Sect. 3.1). 

Given the depth of the imaging data, incompleteness 
is not an issue since our lower magnitude for GC selec- 
tion is ~ 2 mag above the designed photometric limit 
of the images for each galaxy (see Sects. 2.1 and 2.2). 
Further, all dlrr galaxies in our sample have very small 
radial surface brightness gradients which do not bias the 
completeness towards the galaxy center and therefore no 
radial completeness tests are required. 

Although the supreme ACS resolution and deep ex- 
posures allowed us to reject foreground stars and back- 
ground galaxies at high confidence, compact background 
galaxies at intermediate redshifts may share the same 
color, magnitude and size properties as the GCs (e.g. 
Puzia et al. 2004). In order to assess this type of contami- 
nation we used the Hubble Ultra Deep Field (HUDF) im- 
ages. An ideal way to estimate the contamination from 
the HUDF is to run the same image reduction, selec- 
tion and photometry procedures. However, the lack of 
F814W imaging in HUDF prevents us from performing 
this task. Therefore, we selected objects having the same 
characteristics as the GCCs in our F606W images, from 
the HUDF catalog 9 , i.e. 19 < Vp606W < 26 mag, ellip- 
ticity smaller than 0.15 and 2 <FWHM< 9pix. Note 
that Vf606W ~ 26 mag is equal to the magnitude limit 
we used for the most distant galaxies in our sample (i.e. 
UGC3974, UGC3755 at - 8Mpc), therefore the above 
set of parameters for selecting objects from the HUDF 
will be exposure independent. Thus, the objects in the 
HUDF and our fields should be similarly resolved up to 
the magnitude limit used for selection. 

After selecting HUDF objects this way, we converted 
their magnitudes to the Johnson/Cousins system using 
the Sirianni ct al. (2005) transformation coefficients and 
show their magnitude and color distributions in Figure 5. 
The dashed-dotted line, in Fig. 5 a), shows the apparent 
magnitude distribution of all GC candidates in our sam- 
ple while the solid (blue) and dashed (red) histograms 
show the candidates divided by their color into blue and 

9 http:/ /heasarc. gsfc.nasa.gov/W3Browse/all/hubbleudf.html 
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Fig. 6.— 120 X 120 pix F814W excerpts for all GC candidates. 
Labels on the top of the images indicate the name of the host galaxy 
and the sequence number of the GC candidate. Panel a) shows the 
GCs with 0.7 <{V - I)a < 1.15 mag. The last eight GCs, having 
their labels in the bottom of the image, are the faintest and red- 
dest candidates among the blue population (compare with Fig. 7 
and Fig. 9), thus are likely background contaminants. Panel b) 
shows the red ((V — J)o ^ 1.15mag) group of objects. Although 
excluded from the analysis, just for comparison we show U 1281- 
01 and U 3974-05, the brightest objects in this group, having ex- 
tended halos perhaps suggesting that they are background galax- 
ies. This demonstrates that even with the supreme ACS resolution 
only spectral analysis can resolve all doubts regarding their nature, 
especially for objects like U 3974-05. 



red GCCs, respectively (see Sect. 3.1). Wc ovcrplot the 
HUDF objects with a shaded histogram. While covering 
the same (V — I) color distribution as the GC candidates 
(Fig. 5 b) , the expected HUDF contamination from back- 
ground galaxies starts to be a significant factor at faint 
apparent magnitudes, Vq > 24 mag, as seen in Fig. 5 a). 
Hence, we might expect contamination only in the last 
one or two magnitude bins. Given that in Fig. 5 we com- 
bine the GC candidates from 15 individual ACS fields wc 
need to multiply the number of HUDF contaminating ob- 
jects by 15. Considering the apparent magnitude distri- 
bution of the blue and red GC candidates and the prop- 
erly scaled HUDF contaminants, our upper estimates are 
6 and 9 contaminants, respectively. 

In Figure 6 we show 120 x 120 pix F814W image ex- 
cerpts for all GC candidates. With blue and red labels 
in Fig. 6 a), b) are indicated the GCs bluer or redder than 
(V — I)o = 1.1 mag, respectively The objects with la- 
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Fig. 7. — Color-magnitude diagram for all objects that passed the 
GC selection criteria (Sect. 2.2). The colors of the bulk of the GC 
candidates are consistent with old metal-poor GCs. The labels 
indicate the host galaxy and the number of the GC candidates 
in that galaxy. The GC candidates colors and magnitudes were 
dereddened for foreground Galactic extinction using the Schlegel 
ct al. (1998) dust maps. With thick solid and dashed vertical lines 
are indicated the colors separating blue and red GCs. 

bels at the bottom (in magenta) in Fig. 6 a) and b) are 
the reddest GCs with (V - I)o > 1.4 mag (cf. Fig. 7 and 
9) which are very likely background contaminants. Fig. 6 
shows that the probable contaminants are indistinguish- 
able in appearance from the blue GC candidates in the 
ACS images. 

Red GCs are unlikely to be observed in metal-poor dlrr 
galaxies because such red metal-rich GCs are believed to 
be forming mainly in massive galaxies. This question 
will be further addressed in Sect. 3. Therefore, in the 
following sections our analysis will be mainly focused on 
the blue GCCs. 

3. ANALYSIS 

3.1. Magnitudes and Colors 

Given the low luminosities (My > — 17mag, see Ta- 
ble 1) of the dwarf galaxies in our sample and the ex- 
pected low GC specific frequency for dlrr galaxies in 
group environments (e.g. the Local Group, Harris 1991) 
one would expect to find one old GC in each of these 
galaxies. However, our analysis, shows that several of 
the faint dwarf galaxies actually harbor more than one 
GC candidate. As expected, most of the candidates show 
blue (V — I) colors consistent with old metal-poor GCs. 
Figure 7 shows the combined color-magnitude diagram 
for all GC candidates selected among the 19 dlrrs in our 
sample (Sect. 2.2), that fulfill the color, magnitude, size 
criteria and passed the additional visual inspection which 
excluded background galaxies or blends from the final 
list. Other studies (e.g. Seth et al. 2004; Sharina et al. 
2005; Georgiev et al. 2006) of old GCs in dwarf irreg- 
ular galaxies outside the Local Group showed that the 
GCs are typically constrained within the color regime 
0.8 < (V — I)o < 1.1, which corresponds to > 5 Gyr 
old, metal-poor (Z < 0.02 x Zq) GCs according to the 
SSP models (e.g. Bruzual & Chariot 2003). Compared 
with the SSP models, the two very blue {{V - I) < 0.7) 
GCCs (NGC 784-05,06) are likely intermediate-age star 
clusters with an upper age limit of < 2 Gyr. These two 
were therefore excluded from our analysis since we are in- 
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Fig. 8. — Color vs. projected distance for the blue GCs. The iso- 
lines of the two dimensional density estimate shows that there is 
no hint for GCs colors becoming redder with decreasing projected 
distance, as it might be expected due to internal reddening within 
the host galaxy. 

tercstcd only in the old GC population. The population 
of red objects ((V— I)o > 1.1) in our CMD, which passed 
our initial selection criteria, most likely is explained by 
unresolved, redshifted compact galaxies. Further, such 
red, possibly metal-rich GCs, are not expected in dlrrs 
although few of them were observed by Olsen ct al. (2004) 
in Sculptor group late- type dwarfs. It is possible that 
some of the GCs might be reddened, in which case a 
dependence of their colors with the galactocentric dis- 
tance would be expected. It is seen from Figure 8 that 
such a dependence is not present. In fact, we see some 
evidence for a redder color with increasing projected dis- 
tance. Given the low number statistics, this trend could 
be due to unresolved background contaminants in the 
galaxies' outer regions. 

However, follow-up spectroscopy is needed in order 
to constrain their nature. The two bright and red 
objects (My > -7 mag, 1.2 < (V - J) < 1.4) in 
Fig. 7 (UGC 1281-01, UGC 3974-05) show slightly dif- 
ferent (more complex) profiles compared to the rest 
of the GC candidates and therefore are most proba- 
bly background contaminants. We therefore also ex- 
clude these from subsequent analysis. Objects redder 
than (V — I)o > 1.5, although compact and morpholog- 
ically similar to the bluer GC candidates, are deemed 
background contaminants. The structural parameters of 
these objects are discussed in Sect. 3.3. 

A detailed examination of Fig. 7 reveals that there 
seems to be a lack of faint blue GCs in our CMD 
((V — I)o < l;My < —6 mag). We have extensively 
tested the reliability of this observation, by loosening 
our size/shape selection criteria (the limiting IMEXAM 
ellipticity and FWHM f60 6W ~ FWHMpsuw cutoffs). 
This resulted in an increased number of contaminating 
blended stars in crowded regions. Therefore this lack of 
faint bGCs is not due to incompleteness or biased GC 
selection introduced by our size discrimination criteria 
(see Sect. 2.2). However, we do find GC candidates red- 
der (V — I)o ~ 1 and fainter than My ~ —6 mag, which 
are likely background contaminants. This issue is further 
discussed in Sect. 4.1. 

In Figure 9 we present the combined luminosity and 
color distributions of all GC candidates, without the ex- 
cluded objects as discussed above. In order to probe the 
underlying distributions we run a non-parametric prob- 
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Fig. 9. — Color and luminosity functions (panels a and b) of all GCCs detected in 15 dlrr galaxies. The blue histogram marks the 
likely GCs, while the likely background contaminants are plotted as red histogram. With thick and dashed lines are shown the non- 
parametric density estimates using Gaussian and Epanecknikov kernels. Panel b) shows the increasing number of objects in the red group 
with decreasing magnitude, which is mainly due to background contamination (see Sect. 2.7 and Fig. 5). To compare the LF of the GCs in 
our sample with that of the Magellanic-type dlrrs, panels c) and d) show the LF of old GCs in the LMC, SMC and Fornax dwarf galaxy 
(data from McLaughlin & van der Marel 2005). 



ability density estimator using Gaussian and Epanech- 
nikov (inverted parabola) kernel, shown with solid and 
dashed lines in Fig. 9 a). Although the Gaussian kernel 
shows better the dip at (V — I)q ~ 1.1 between the red 
and blue objects, both kernels give similar color number 
density probability distributions, with maxima and min- 
ima at the same colors. The third peak represents the 
population of very red objects, i.e. contaminating back- 
ground galaxies. The highest probability density peak 
values are at (V ~ I) a ~ 0.96 and (V - 7) ^ 1.19 for 
the blue (bGC) and red GCs, respectively. The peak 
(V — I) color of the bGCs is similar to the one found in 
other dlrr galaxies (e.g. Seth et al. 2004; Sharina et al. 
2005; Georgiev et al. 2006). Although our GC sample 
might have reddening uncertainties, it is anticipated that 
it is not a significant factor for these extremely low-mass 
galaxies. Based on H a fluxes James et al. (2005) showed 
that dlrrs posses low Au a — 0.5, which translates to 
Eb-v = 0.06 or Ey-i = 0.07 mag. Note, however, that 
this only refers to the star forming regions. 

The luminosity distribution of the bGCs in our sam- 
ple peaks at Mto,v = —7.43 (rather typical for mas- 
sive galaxies) and has a broader dispersion than the ob- 
served LF for the LMC, SMC and Fornax dwarf galaxies 
(data from McLaughlin & van der Marel 2005), shown 



in Fig. 9 c) . The broader luminosity distributions might 
be partly attributed due to the distance uncertainties to 
these dwarfs of a(m - M) ~ 0.05 (Tully et al. 2006). A 
Kolmogorov-Smirnov test gives a 49% probability that 
the My distribution of our bGCs are drawn from those of 
the LMC. In comparison with other low-mass dlrr galax- 
ies (Mto,v ~ —6; Sharina et al. 2005; van den Bergh 
2006) the bGCs in our sample show a brighter Mto- 

It is interesting to note that the LMC data from 
McLaughlin & van der Marel (2005), which is the re- 
calibrated and extended version of the LMC data from 
Mackey & Gilmore (2003a, b), gives a brighter peak mag- 
nitude at Mto.v = — 7.50 ± 0.16 mag (highest probabil- 
ity density peak at My = —7.52 mag) versus Mto.v = 
—7.31 mag, derived from fitting the formerly determined 
M v values of old LMC GCs (Table 2 in van den Bergh 
& Mackey 2004). Given the uncertainties, this suggests 
that the GCLF of old GCs for the dlrrs in our sample, 
and LMC (considering the uncertainties), both peak at 
similar magnitudes. 

Much fainter Mto values were found for other dlrrs 
by Sharina et al. (2005) using WFPC2 data. However, 
some of the galaxies included in their study (eg. KK 16, 
UGC8651, UGC4115, UGC3974, UGC3755) were re- 
observed with much deeper ACS observations by Tully 



10 



Gcorgiev, I. Y. et al. 





-0.1 




-0.2 




-0.3 




-0.4 


o 

> 


-0.5 


AM 


-0.6 




-0.7 




-0.8 




-0.9 




-1 



> 
< 



a) 



(This work) - (Sharina'05) 



U37f5-01 U37f5-04 



U37f5-02 U37!5-( 



AM V0 =-0.52 +0.03 



-9.5 



8.5 -8 -7.5 
My (This work) 



-6.5 



0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 



-0.5 



b) 



U37:5TO2 
IB 7i 5-04 



A(V-I) =-0.07 +0.04 



0.75 0.8 0.85 0.9 0.95 1 
(V-I) (This work) 



1.05 1.1 



Total (150) 
YH (30) 
OH (70) 
BD(37) ■ 

= -7.27 ±0.11 




-11 -10 -9 



Total (95) - 
YH(17) KJ 
OH (61) 
BD(17) 



0.6 0.7 0.8 0.9 




Fig. 10. — Differences in absolute magnitudes and colors for 
our GCs in common with Sharina et al. (2005). The least-square 
fit to the data (solid line) shows that our magnitudes are 0.5 mag 
brighter than the earlier study, while the expected difference due 
to distance modulus is -0.76 mag (dotted line). The colors (within 
the errors) are basically identical between the two studies. 

et al. (2006). For these dlrr the distance moduli were im- 
proved and some of them turned out to be more distant 
than thought before. Therefore, underestimating the dis- 
tance to some galaxies of the Sharina et al. (2005) sample 
led to an underestimate of the GCs' absolute magnitudes, 
and hence the Mto magnitude. The larger distances also 
have an effect on the GC half-light radius estimates (see 
Sect. 3.3). For the GCs in common with the Sharina 
et al. (2005) study in UGC 3755, the differences in abso- 
lute magnitude between both studies are mainly due to 
the improved photometric technique we used for deriving 
the clusters apparent magnitudes and the new distance 
moduli. Figure 10 shows the differences between both 
studies. However, the faintness and broadness of the 
Sharina et al. (2005) GC LFs is not expected to be influ- 
enced by contamination from background galaxies since 
their estimated contamination level is less than ~ 5% for 
objects with (V - I) < 1.2, -5.3 < M Vo < -6.3 mag 
(Puzia & Sharina 2008). 

3.2. Comparison with the Milky Way GCs 

In Figure 11 we show the color and magnitude distribu- 
tions of the different GC subpopulations in our Galaxy. 
The Milky Way "Young Halo" (YH) GCs, which were 
suggested to have been accreted from dwarf galaxies (e.g. 



Fig. 11. — Color and magnitude distributions of Galactic GCs 
(data from Harris 1996, revised 2003). Based on their metallicity 
and horizontal branch morphology, Bulge/Disk (BD), Old Halo 
(OH) and "Young" Halo (YH) GCs are shown in red, green and 
blue, respectively. Excluding the last two bins, a Gaussian fit to 
the YH GC histogram gives a Mto YH value that is fainter than 
those of our dlrrs' GCs and the LMC GCs (see Fig. 9). 



Zinn 1993; Mackey & Gilmore 2004; Mackey & van den 
Bcrgh 2005; Lee et al. 2007) show a turnover magni- 
tude {Mto = —7.27 ±0.11 mag) which, although fainter, 
is consistent within the uncertainties with that of our 
bGCs (Mto = —7.41 ± 0.22 mag). In comparison with 
the Magellanic Clouds (see also Sect. 3.1), the magni- 
tude distribution of the LMC old GCs peaks at surpris- 
ingly bright magnitude M T o = -7.78 ± 0.12. When 
including the old GCs of the SMC and the Fornax dSph 
the combined peak magnitude becomes slightly fainter 
(Mto = -7.67±0.1mag). K-S tests (cf. Table4) assign 
54% and 50% similarity for the V — I and My distribu- 
tions, respectively, between our bGCs and the Galactic 
YH clusters. It should be noted, that only 17 YH GCs 
have V — I colors in the Harris (1996) catalog. Com- 
paring with the same distributions for the OH GCs gives 
probabilities of 2% and 19%. The luminosity function of 
the old LMC clusters is similar to the YH and OH clus- 
ters with K-S probability of 43% and 78%, respectively. 
In Table 4 are also shown the results from K-S tests be- 
tween the V — I and My distributions for ten of the old 
GCs in the LMC with V-I colors (McLaughlin & van der 
Marcl 2005) and bGC in our dlrrs sample. 



Old Globular Clusters in Magellanic- Type dlrrs. 



11 



64 
32 
16 
8 
4 
2 
1 



a ) bGCs 
CO Cen 
Red 

log(r h ) = 0.25xM v + 2.95 



-11 



_jbf959-04 



- " *nw(m-03 W $mm V 959 ' 02 # 



-10 



-7 



-6 



14 
12 
10 

8 
6 
4 
2 





bGCs (37) rxWxM 



12 15 18 21 24 



11 

9 
7 
5 
3 
1 



-b) 

784-04 


bGCs ^ . 
CoCen 

Red ' ® ' _ 


^1311-06 | 

THil54-023-0] 


t 


f *: 1959-07 

fit 1 V V"^ 9sw " 


AtK 16-01 , 
fJLMl-04j| 3 ii_oi 

$1959-06 § 



0.5 



1.5 



2.5 



3.5 



d prq/ r eff 



LMC + SMC + Fnx (22) 

b) LMC (16) i xxxx> < 




12 15 18 21 24 



Fig. 12. — a) Distribution of the GC candidate sizes (r^) with 
their absolute magnitude, b) rj, as function of the normalized 
projected distance (d pro j /V e fl ) from the host galaxy center (b). 
The GC candidates occupy the My vs. plane below the relation 
defined by Mackey &: van den Bergh (2005) (thick line in panel a), 
which is typical for GCs in our Galaxy and other galaxies. IC 1959- 
04, the nuclear cluster of IC 1959, shares the same parameter space 
as ijj Cen. 

3.3. Structural parameters 
3.3.1. Half-light radii 

The half- light radius r^, a measure of the GC size and 
stable over many (> 10) relaxations times (e.g. Spitzer & 
Thuan 1972; Aarseth & Heggie 1998), is the most robust 
and easy to measure structural parameter and therefore 
is a good indicator for the initial conditions of the GC 
formation. Hence, a comparison between r\ L of the blue 
GCs in our dlrr sample and r^ of the YH GCs in our 
Galaxy and in the LMC can reveal an important evolu- 
tionary link among them as well as information of the 
initial conditions of GC formation. 

The rh versus My relation in Fig. 12 a) shows that most 
of our GC candidates fall below the upper envelope of 
the GC distribution in our Galaxy (solid line in Fig. 12; 
see Mackey & van den Bergh 2005). Interestingly, the 
brightest cluster, IC1959-04, is at the very center of its 
host galaxy and occupies the same r^ versus My region 
and has similar color as to Cen, which supports a sce- 
nario for u> Cen's origin as the nucleus of a former dwarf 
galaxy (e.g. Hilker & Richtler 2000; Hughes & Waller- 
stein 2000). However, this is the only nuclear cluster we 
find in the studied galaxies. It is known that r^ increases 
with galactocentric distance for MW, M31 and Cen A 
GCs (e.g. Barmby et al. 2007), which indicates the influ- 
ence of the galactic potential. In Fig. 12 b) we combined 
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Fig. 13. — Comparison of ri, distributions for (a) GCs in our 
dlrrs with that of (b) LMC (+SMC+Fornax) GCs and (c) MW 
GCs (data from Mackey & van den Bergh 2005; McLaughlin & 
van der Marel 2005). 



the GCs rhS from all galaxies in our sample. Note that 
in Fig. 12 b we define r e g = y/ax b, where a and b are 
the galaxy semi-major and semi-minor axes of the galaxy 
light profile. A dependence of rh on the projected dis- 
tance from the galaxy center can not be confirmed with 
significance given the large scatter and the low number 
statistics. 

In Figure 13 we show the size distribution of the GCs 
in our sample compared with that of the old GCs in 
the Magellanic Clouds and GCs in our Galaxy. Our 
bGCs are more compact on average (fh ^3.3 pc, cf. 
Fig. 13) than the LMC GCs (f h - 6pc) and the YH 
GCs (fh ~ 7.7 pc). They are rather more similar to the 
MW OH GCs (f h ~3.5pc), as evidenced by the K-S test 
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GCs in common with Sharina ct al. (2005). The least-square fit to 
the data (solid line) shows that our sizes are ~ 5 pc smaller than 
found in the earlier study, while the expected difference due to a 
different distance modulus is only — 0.55pc. 

probabilities listed in Table 4. In Sect. 4.2 we discuss the 
implication of this comparison. 

It is worth noting, that in the SMC, a dlrr galaxy 
which closely resembles in total luminosity (My,SMC = 
— 16.82 mag) the galaxies in the current study, there ex- 
ists only one old GC (10.6 Gyr, Dolphin ct al. 2001) with 
rh — 5.6 pc. 

Probability K-S tests give 1% and 8% likelihood that 
the bGC have a similar distribution as the YHs or the 
old LMC GCs, respectively. The r h distribution of LMC 
GCs is 32% and 7% similar to the ones of YHs and OHs, 
respectively (see Table 4). 

GCs more extended than those of the Milky Way and 
elliptical galaxies have been previously observed in other 
dlrrs (Seth et al. 2004; Sharina et al. 2005). We have 
nine GCs in common (one in UGC4115 and eight in 
UGC 3755) with the Sharina et al. (2005) study. Our r h 
measurements, which take into account the cluster ellip- 
ticity, are significantly smaller by ~ 5 pc (see Fig. 14) . 
The expected difference due to updated distances to 
UGC 3755 (m-M= 29.35, Tully et al. 2006) between 
both studies is just — 0.53pc (m-M= 28.59, Karachcnt- 
sev et al. 2004, adopted by Sharina et al. (2005)). 

However, the absolute rh comparisons between GCs 
should be taken with caution as they were assessed using 



different measurement methods by different studies. The 
structural parameters of the Galactic and LMC GCs are 
based on King model fits from McLaughlin & van der 
Marcl (2005). The same King model was used in our 
study, but using discrete concentration parameters as 
provided by ishape, while McLaughlin & van der Marel 
(2005) derived a concentration parameter coming from 
their best \ 2 fit to the model. 

3.3.2. Ellipticities 

Two major dynamical effects are considered as drivers 
of the evolution of the GC ellipticities (e): i) the external 
galactic tidal field and/or ii) internal (age-dependent) 
rotation and/or velocity anisotropy (e.g. Fall & Frenk 
1985; Han & Ryden 1994; Meylan & Heggie 1997). A 
third possibility, although perhaps less attractive than 
general mechanisms, is that highly flattened GCs might 
arise from cluster mergers (e.g. de Oliveira et al. 2000). 
Fall & Frenk (1985) showed that the cluster e decreases 
by a factor ~2-5x t r h- The relaxation time depends 
on the clusters mass at r = r^. For Galactic GCs it 
ranges from ~ 10 8 to ~ 10 10 yr (Djorgovski 1993) and 

- 10 8 10 9 yr for LMC GCs (Fall & Frenk 1985). It 

is well known, that on average the Magellanic Clouds 
GCs are flatter than those in the Milky Way and M31 
(e.g. Geisler & Hodge 1980) and these differences might 
be primarily due to differences in cluster ages (Han & 
Ryden 1994). 

The ellipticities of the clusters in our sample are mea- 
sured at the rh along the semi-major axis returned by 
ISHAPE. In Figure 15 we present the ellipticities of the 
bGCs in our dlrr galaxies versus My, dp ro j/r ff , (V — I)o 
and in panels a,b,c and d, respectively. For com- 
parison we show the old GCs in the Magellanic Clouds 
with data available from the literature (Frenk & Fall 
1982; Kontizas et al. 1989, 1990). In general, the old 
LMC/SMC GCs overlap well with the clusters in our 
sample. As previous studies for various galaxy types also 
showed, no correlation between e and the projected dis- 
tance from the galactic center can be seen (Fig. 15 b). 

The comparison between the ellipticity vs. My, (V — 
I)o and rh distributions the bGCs in this study with 
those of YH GCs in our Galaxy (suspected of being ac- 
creted from dwarf galaxies), is shown in Figure 16. For 
the Galactic clusters, we used the e entries in the 2003 
update of the Harris (1996) catalog. Although only 16 
YH GCs are included, it is seen that they overlap with 
the old Magellanic GCs in the e vs. My and r\ x planes 
(Fig. 16a, c). However, the e distributions (Fig. 17) of 
our bGCs and the LMC GCs are very similar to one an- 
other (e ~ 0.12 ± 0.07 and e ~ 0.11 ± 0.06) and flatter 
than that of the YH GCs (e ~ 0.06 ± 0.06), and show 
a peaked, Gaussian-type distribution. The result of a 
Kolmogorov-Smirnov test (see Table 4) run on the e dis- 
tributions of the YH GCs, LMC old GCs and the bGCs 
in our sample showed that there is only 1 and 3% chance 
that the LMC old GCs or bGCs, respectively, are drawn 
from the same distribution as the YH GCs, while the 
bGCs and the LMC GCs share the same distribution at 
a 99% confidence level. The presence of such a turnover 
might mean that either the clusters are not fully relaxed, 
i.e. they are dynamically young and still evolving, or 
it is an indication of the impact of the weaker potential 
of these dwarfs on the initial conditions during the GC 
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formation, or both. 

Therefore, if the YH GCs are of external origin, their 
e distribution should be similar to that of the bGCs. In 
contrast, it is more like a power law as seen in the Fig. 17, 
lower panel. This might be an indication that their accre- 
tion must have happened very early and the strong tidal 
field of our Galaxy might be responsible for the evolu- 
tion of the e distribution away from the initial conditions, 
removing the turnover. Note, however, that reliable el- 
lipticities are available only for approximately half of all 
MW YH GCs (16 out of 30) and for old GCs in the LMC 
(11 out of 17). Thus the ellipticity distribution we show 
might be biased. 

3.4. Specific frequencies 

The 'specific frequency' (Sn) is a quantity introduced 
by Harris & van den Bergh (1981) to compare the rich- 
ness of GCSs between (elliptical) galaxies. The Sn was 
originally defined as twice the number of GCs down to 
the turnover peak of the GC luminosity function, or (as 
most commonly used) as the total number of GCs (Nqc), 
normalized to a galaxy luminosity of My = —15 mag: 

S N = N GC 10 OA ( Mv+15 > (2) 

For poor GCSs containing few GCs, such as those in 
our dlrrs, calculating the S n as originally defined is prac- 
tically impossible. Hence, the Sn values presented here 



were derived from the total Nqq- Given the small phys- 
ical sizes of the studied dlrrs and the large enough ACS 
field of view, there is no need to apply a geometrical 
completeness correction to the total GC number for most 
of the galaxies except for the brightest and largest two, 
NGC784 and ESO 154-023. Because of the extremely 
deep ACS observations, corrections due to photometric 
incompleteness arc not necessary either (see Sect. 2.7). 
In order to derive the total galaxy light we used the 
IRAF/STSDAS task ellipse. The galaxy magnitudes 
we derived (see Table 1) were identical (within the errors) 
to the RC3 (de Vaucouleurs et al. 1991) values for the 
same galaxies at the radius of the {1b — 25 mag arcsec -2 
isophote. However, we could not directly measure the 
total magnitudes for NGC 784 and ESO 154-023 which 
extend beyond the ACS field and thus we adopted their 
RC3 magnitudes. 

The majority of the bGCs in our ACS images are re- 
solved and the estimated background contamination for 
this GC sub-population is low (see Sect. 2.7). Thus, the 
uncertainty of the Sn values for the bGCs presented 
in Table 2 is mainly due to stochasticity in the num- 
ber of genuine GCs rather than background contami- 
nation. The subpopulation of red GCs is more likely 
affected by faint and unresolved background sources, be- 
cause red GCs are on average fainter and more compact 
than bGCs. This is the reason why our analysis is mainly 
focused on bGCs. 
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ues of the Magellanic Clouds and NGC 1427A are shown 
for comparison. To calculate the Magellanic Clouds' ab- 
solute magnitudes we used the distance moduli of Lancy 
& Stobie (1994) and information for NGC 1427A was 
taken from Georgiev et al. (2006). Note that the dlrrs 
in our sample with luminosities similar to the SMC also 
have similar Sn values. In general we find that the GC 
specific frequency in dlrrs increases (on average) with 
decreasing host galaxy luminosity. This result was also 
found for early- type dwarfs in galaxy clusters (Miller 
et al. 1998; Forbes 2005; Strader et al. 2006; Miller & 
Lotz 2007). 



3.4.1. Present-day Sn Values 

Using the distance modulus for each galaxy as deter- 
mined by Tully et al. (2006) the apparent galaxy mag- 
nitudes were converted to absolute ones. Taking the 
numbers of GCs as given in Table 2, the resulting Sn 
values then range from 0.74 to 88. Excluding the ex- 
treme outlier (KK 27) with Sn = 88, the mean value for 
all galaxies with blue GCs is ~ 3.8. The total mean Sn 
would increase to ~ 4.2 if one would include the red GCs. 
Figure 18 shows the Sn values for each individual galaxy 
(symbols) and an average Sn for galaxies in magnitude 
bins (histogram). The galaxies in which no GC was de- 
tected are plotted with negative Sn values. The Sn val- 



3.4.2. Sn Values After Evolutionary Fading 

Since the dlrr galaxies in our sample are of Magellanic 
type with ongoing star formation, their present-day in- 
tegrated luminosity is dominated by the young stellar 
populations that formed in the most recent star forma- 
tion episode. In order to make a more relevant compar- 
ison with the Sn values in evolved dwarf galaxies like 
dEs or dSphs, we have to account for fading of the dlrrs' 
luminosity due to stellar evolution. All galaxies in our 
sample have luminosities comparable to or even fainter 
than the SMC. We therefore adopted the SMC metallic- 
ity of [Fe/H]= -0.68 (Luck et al. 1998) for the passive 
evolution of the stellar population. Bruzual & Chariot 
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(2003) SSP models were used to derive the stellar age for 
each galaxy based on their color and adopted metallicity. 
The mean stellar age for all galaxies as inferred by their 
V — I colors is 1.3 Gyr, with individual ages ranging from 
0.4 to 3.7 Gyr. Given the low surface brightness and low 
star formation of these dwarfs, a passive (dissipationless) 
evolution was assumed to evolve the galaxies' luminosity 
to an age of 14 Gyr. This resulted in an average luminos- 
ity fading of 2.34 ± 0.43 mag, with a full range from 2.98 
to 1.09 mag. This leads to an increase in the Sn values 
by a factor of 2 to almost 16 due to passive evolutionary 
fading. The specific frequencies after age fading (Sn/) 
are listed in Table 2. 

It is worth assessing the impact of the adopted SMC 
metallicity on the evolved Sn values. Due to the contin- 
uous star formation rate, typical for late-type galaxies, 
a dlrr galaxy in our sample cannot be more metal-poor 
than its most metal-poor GC(s); unless, in the unlikely 
scenario where the galaxy underwent instantaneous star 
formation which ceased abruptly prior to GC formation. 
However, assuming a mean galaxy metallicity as metal- 
poor as 0.02 x Zq, the mean stellar age, as inferred from 
the measured V — I color, increases by ~ 1 Gyr. Hence, 
after evolutionary age fading the galaxy will appear, on 
average, ~ 1 mag fainter. This translates into an increase 
of Sn up to 6 times of the present-day value, unlike up 
to 16 times if the dlrrs were of SMC metallicity. In gen- 
eral, the change of Sn values due to evolutionary age 
fading will be milder for more metal-poor galaxies, and 
the evolved S n values in Table 2 will apply only to galax- 
ies with SMC metallicity ([Fe/H] = —0.68, i.e. similar to 
0.2 x Z Q in the Bruzual & Chariot (2003) SSPs). This 
illustrates the dependence of Sn values on host galaxy 



metallicity. 

If one would add all galaxies with GCs together, a 
master galaxy of My = —18.24 mag would result with 
a SNb = 1-9. Including also dlrrs without GCs results 
in a master galaxy of My = —18.36 mag and SNb = 1-7. 
Fading the master galaxy light to an age of 14 Gyr gives 
for the two cases My = —16.0 and My = —16.11 mag 
and specific frequencies of Sn = 14.8 and Sn = 13, re- 
spectively. Again, this is consistent with Sn values for 
dEs/dSphs of similar absolute magnitude. 

4. DISCUSSION AND CONCLUSIONS 

We have searched for old globular clusters in 19 
Magellanic-type dwarf Irregular (dlrr) galaxies utilizing 
archival F606W and F8UW HST/ACS images. The 
main goal of this study is to assess the properties of the 
old GCs in such systems and compare them with those 
of the Galactic "Young Halo" GCs (YH GCs), which 
are considered to be of external origin (e.g. Zinn 1993; 
Mackey & van den Bergh 2005; Lee et al. 2007). The 
studied dlrrs reside in nearby (2-8 Mpc) groups and as- 
sociations containing dwarf galaxies only, with no mas- 
sive nearby galaxy around (Tully et al. 2006). All galax- 
ies in our sample have absolute magnitudes smaller than 
or equal to the SMC. 

4.1. An Apparent Lack of Low-Luminosity GCs 

We detect in total 50 GCs in 13 dlrrs, of which 37 have 
(V — I) colors consistent with "blue" (old and metal- 
poor) GCs (bGC) typical for low-mass dwarfs and the 
bGC population in other galaxies. The few red GCs 
(rGC) we detect are indistinguishable in sizes from the 
bGCs and are not expected in this galaxy type, perhaps 
they are background contaminants. Only follow-up spec- 
troscopy can reveal their nature. The luminosity func- 
tion (LF) of the bGC candidates in our sample shows 
a turnover magnitude at My = —7.41 ± 0.22 mag, simi- 
lar to what is typically found in other galaxy types, and 
a GCLF width of a = 1.79 ± 0.31, typical for dlrrs but 
slightly broader than the typical width of GCLFs in mas- 
sive galaxies. The broadening of the GCLF might be (at 
least partly) due to the distance modulus uncertainties 
to those dwarfs of cr(m - M) ~ 0.05 (Tully et al. 2006). 

Interestingly, we do not find many low-luminosity GCs, 
contrary to expectations. As we showed, this is not due 
to incompleteness or selection criteria (Sect. 3.1). The 
lack of faint GCs might be related to environmental ef- 
fects and the physical conditions of the GC formation 
and disruption in such dwarfs. However, all dlrrs in our 
sample are found in associations of dwarf galaxies, e.g. 
there is no dominant massive galaxy in their immediate 
vicinity. Thus tidal stripping or disruption in the current 
environment is unlikely to be responsible. 

Another possibility could be that the bluer (0.8 ~ 
V — I) and slightly brighter (My ~ —6.5 mag) GCs seen 
in Fig. 7 might actually be relatively young GCs. Aging 
due to stellar evolution would make their magnitudes 
fainter and colors redder. If so, a shift in their colors 
and magnitudes of (V — I) ~ 0.15 and My ~ 1.5 mag 
would be required to populate that region in the CMD. 
According to the Bruzual & Chariot (2003) SSP models, 
this is exactly the amount of color and magnitude evo- 
lution needed for a > 4 Gyr old metal-poor cluster with 
(V - I) > 0.8, to evolve to ~ 14 Gyr. 
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Further detailed spectral analysis of GC candidates in 
dlrrs are necessary to address this question and draw 
firmer conclusions. Although an interesting observa- 
tional result, for now we can not extend the discussion 
further given the inadequate number statistics in our 
data set. 

4.2. Sizes: GCs in dlrr Galaxies are smaller than 
"Young Halo" GCs 

The half-light radii, rh, of the bGCs in our sample are 
smaller on average, rh ~ 3.3 pc, than those of the YH 
GCs in our Galaxy (rh ~ 7.7 pc) and those in the Magel- 
lanic Clouds and Fornax dwarfs (rh ~ 6.3 pc) and other 
dlrrs (Sharina et al. 2005; van den Bergh 2006). They 
compare well with the rh values of Galactic OH GCs, as 
shown by the K-S tests in Table 4. This suggests that 
either the conditions of YH GC formation were differ- 
ent from those in low-mass galaxies, or their sizes were 
influenced by tidal shocks in their subsequent evolution 
in the Galactic tidal field which will require that they 
have been accreted very early. N— body simulations (e.g. 
Dinescu et al. 1999; Gieles et al. 2007) show that the 
effect of bulge/disk shocks scales with r^/M c \. Hence, 
depending on their orbital periods and pericentric radii, 
the tidal field will strongly affect diffuse clusters, such as 
the YH GCs, and transform their rh distributions. Fur- 
ther, Lamcrs et al. (2005) show that the cluster disrup- 
tion time scales with the ambient density. As rh is stable 
over many relaxation times (Sect. 3.3.1) in the absence 
of a tidal field, a modeling of the rh evolution with the 
strength of the host galaxy tidal field and environmental 
density will be required to understand the rh evolution. 

4.3. The Ellipticity Distribution of GCs in dlrr 
Galaxies 

We find that, on average, the ellipticities (e) for the 
GCs in the studied dlrrs are practically the same as those 
of the old GCs in the Magellanic Clouds (e ~ 0.12 ± 0.07 
and e~ 0.11±0.06, respectively). They also show sim- 
ilar distributions of e with My, rh and cluster color 
(Fig. 15). Although ellipticities are only available for half 
the Galactic YH GCs, their distribution differs signifi- 
cantly from those of GCs in our sample and those in the 
LMC (see Fig. 17). This is further shown by the result 
of the K-S tests presented in Table 4. This shows that if 
the YH GCs formed in similar environments as the GCs 
in low-mass dlrr galaxies, their ellipticity must have been 
significantly altered by the tidal field of the Galaxy after 
their accretion. I.e., if YH GCs were accreted from satel- 
lite galaxies, the accretion event (s) must have happened 
quite early 

The observed turnover in the ellipticity distributions 
for the bGCs and LMC old GCs in contrast with the 
power-law e distribution for the Galactic YHs, might sug- 
gest that the former are slightly younger and not yet fully 
relaxed (see also Sect. 3.3.2), as expected from theory 
(Fall & Frenk 1985). 

Furthermore, the strength of the galactic tidal field 
is expected to be important in setting the initial condi- 
tions during the GC formation. It probably also causes 
the mean e value to vary with galaxy mass, as slightly 
higher e values are observed for the young clusters in the 
SMC than in the LMC (Kontizas et al. 1990). Thus, it 
does not seem surprising that the GCs in our low-mass 



dlrrs and in the Magellanic Clouds on average are flat- 
ter than those in massive galaxies and that, if the YH 
GCs were accreted early enough, their evolution in the 
stronger Galactic tidal field must have influenced their 
ellipticities. The fact that the evs.rh relation for the 
YH GCs holds just as well for GCs in our dlrrs and the 
LMC/SMC suggests that they might all have formed in 
similar initial conditions. However, this interpretation 
needs to be tested with larger numbers of GCs in dlrrs, 
simulations and further precise determinations of the GC 
ellipticities in various environments, as well as future in- 
tegrated spectra to infer their ages and metallicities. 

4.4. Specific Frequencies: Trend with galaxy luminosity 
and Comparison to Dwarf Ellipticals 

The present-day specific frequencies (Sn) of the dlrr 
galaxies in our sample span a broad range from Sn = 0.3 
to 87. Using the galaxy V — I colors and the SSP models 
of Bruzual & Chariot (2003) with the metallicity of the 
SMC, we derive age fading of the galaxies' luminosity to 
14Gyr. We find that this causes an increase in the Sn 
by a factor of ~ 2 to 16, if all GCs survive and no new 
ones are formed (see also Table 2). 

Assuming that these groups of dlrr galaxies might not 
be fully relaxed yet (Tully et al. 2006) and that their 
members will merge eventually, we constructed a "mas- 
ter" dlrr galaxy from all members. The resulting present- 
day and faded Sn values of the "master" dlrr are similar 
to the values for early-type dwarfs. This supports the 
idea that some dE/dSphs might have evolved from dlrrs 
(Miller et al. 1998), due to gas loss after the main burst 
of star formation (Dekel & Silk 1986) and/or via external 
ram pressure stripping (Moore et al. 1998; Grebel et al. 
2003). The Sn values of the non-faded and faded dlrrs, 
in general, are consistent with those found for dEs in the 
same luminosity range (Miller et al. 1998; Miller & Lotz 
2007). 

The recent discovery of hidden non-axisymmctric disk- 
like structures and other morphological anomalies in 
Virgo early-type dwarfs from detailed isophotal analysis 
(Ferrarese et al. 2006) also suggests that their progenitors 
might be late-type dwarfs. This view is theoretically sup- 
ported by Mastropietro et al. (2005). If rapid cessation of 
star formation processes occurs in present-day dlrrs, the 
comparison with Bruzual & Chariot (2003) stellar evolu- 
tion models shows that ~ 10 Gyr of passive fading are re- 
quired for dlrrs to fade and fall on the dSphs metallicity- 
luminosity relation (see also Jerjen et al. 2004). 

If all dlrrs in our sample are merged together (without 
dissipation) to a master dlrr galaxy and the galaxy light 
is evolved to 14 Gyr, the GCS characteristics (Table 2) 
of the "merger remnant" resemble those of a field dwarf 
elliptical or spheroidal galaxy as suggested by Miller et al. 
(1998). 

5. SUMMARY AND OUTLOOK 

Our study shows that properties of old GCs in low- 
mass dlrr galaxies hold interesting and important in- 
sights in our understanding of what were the conditions 
at the time of their formation. Indeed, considering their 
color and magnitude distributions, they are as old and as 
metal poor as the GCs in other galaxy types. Their struc- 
tural parameters (r^ and e) , whose evolution is shown to 
depend mainly on processes internal to the GCs (e.g. ro- 
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tation, two-body relaxation), are very similar to those of 
the old GCs in the Magellanic Clouds, but differ from YH 
GCs in our Galaxy. Under the assumption that initial 
conditions of GC formation (occurring in giant molecu- 
lar clouds) were similar among old GCs, this difference 
suggests that the initial structural parameters were al- 
tered by an amount depending on the environment, i.e., 
the strength of the galactic potential plus external pres- 
sure. Therefore, our comparison of the bGCs in dlrrs 
with the Galactic YH GCs, the latter of which are sus- 
pected to have their origin in accreted low-mass galaxies, 
may indicate that the Galactic YH GCs were accreted 
very early. This accretion must have been early enough 
for the Galactic tidal field to have an impact on the time 
scale of the evolution of the structural parameters. 

However, to properly assess all the addressed sugges- 
tions to a higher statistical confidence, we suggest that 
further detailed studies of GC properties in other sam- 
ples of low-mass dlrr galaxies be undertaken in the near 



future. 
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TABLE 1 

General properties of the studied dwarf irregular galaxies. In column (1) are listed the Group ID to which the galaxy 
belongs and the galaxy's name; columns (2) and (3) list the galaxies' coordinates, columns (4) and (5) their morphological 

classification; columns (6) and (7) give the distance and distance modulus from tully et al. (2006), followed by the 
galaxies' apparent (8) and absolute (9) magnitudes and colors (10) measured in this work (all dereddened from galactic 

extinction). 
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a Tully et al. (2006) 

b Information for UGC 9240, which belongs to this group, was taken from Sharina ct al. (2005) 
c From RC3 
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TABLE 2 

Specific frequencies. In column (1) are listed the Group ID to which the galaxy belongs and the galaxy's 

NAME; THE FOLLOWING COLUMNS ARE: THE DISTANCE (2), ABSOLUTE MAGNITUDE (3), NUMBER OF BLUE GCS (4), NUMBER 

OF RED GCCS (5), BLUE GCS AND TOTAL SPECIFIC FREQUENCY Sjv;, (5jvt D t) (6), AGE OF THE DOMINANT STELLAR 
POPULATION DEDUCED FROM THE BRUZUAL &; CHARLOT (2003) SSP (0.2 X Zq) USING THE GALAXY COLOR (7), AMOUNT OF 
V — BAND FADING TO 14 GYR (8), Myj AFTER FADING (9), AFTER FADING (10), My OF THE PRODUCT GALAXY AFTER 

MERGING ALL GROUP MEMBERS TOGETHER (11), SjV OF THE MERGER PRODUCT, AMOUNT OF THE MERGER PRODUCT FADING 
TO 14 GYR. (12), M VgIp j AFTER MERGER FADING (13), Sfj AFTER MERGER FADING (14). NUMBERS IN BRACKETS ALWAYS 

REFER TO THE TOTAL VALUES. 
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7.41 


-15.50 


8 


1 


5.05 


1.3 


2.43 


-13.07 


47.3 






















(5.68) 








(53.3) 












Dregs 






























UGC 9128 


2.24 


-12.06 








(0) 


0.5 


2.85 


-9.21 


(0) 













a Tully et al. (2006) 

b Information for UGC 9240, which belongs to this group, was taken from Sharina et al. (2005) 
c From RC3 
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TABLE 3 

General properties of the GC candidates in the studied dIrrs. According to the colors of the GC 

CANDIDATES, THEY ARE SEPARATED WITH LINES ON BLUE, RED AND LIKELY BACKGROUND CONTAMINATION SUB-POPULATIONS. 



ID X , Y R.A DEC (2000) MZ (V - Z)o r h r h e d pro j d£ 





(P 


x) 


(hh:mm:ss) 


(mag) 






(mag) 




(pix) 






(pc) 






(kpc) 


(kpc) 


E154-023-01 


3472.05 


2329.00 


02:56:49.17 


-54:33:16.51 


-4.78 ± 


04 


1 


05 ± 0.01 


2 


02 ± 


09 


4 


68 ± 0.15 





10 


1 


99 


0.75 


IC1959-02 


3083.73 


1739.38 


03:33:14.24 


-50:25:40.16 


-6.19 ± 


05 


1 


03 ± 0.01 


4 


23 ± 


01 


6 


21 ± 0.11 





29 


1 


71 


1.42 


IC1959-04 


2224.47 


2273.96 


03:33:12.43 


-50:24:52.60 


-9.83 ± 


05 





97 ± 0.01 


2 


11 ± 


00 


3 


27 ± 0.10 





OS 





23 


0.19 


IC1959-05 


1982.16 


2506.51 


03:33:11.51 


-50:24:38.34 


-6.54 ± 


05 


1 


10 ± 0.01 


1 


81 ± 


00 


2 


81 ± 0.10 





15 





20 


0.22 


IC1959-06 


544.56 , 


3554.56 


03:33:07.70 


-50:23:17.16 


-7.60 ± 


05 





90 ± 0.01 


1 


72 ± 


02 


2 


67 ± 0.12 





11 


2 


87 


2.39 


IC1959-07 


2852.76 


1501.31 


03:33:15.71 


-50:25:31.32 


-6.11 ± 


06 





78 ± 0.03 


2 


66 ± 


05 


4 


12 ± 0.15 





06 


1 


60 


1.41 


IC1959-08 


3366.00 


2605.55 


03:33:09.50 


-50:25:45.00 


-4.65 ± 


05 





99 ± 0.01 


1 


40 ± 


07 


2 


17 ± 0.17 





04 


1 


87 


1.56 


IC1959-09 


2884.04 


3886.73 


03:33:03.47 


-50:25:08.15 


—5.06 ± 


06 


1 


08 ± 0.02 


1 


67 ± 


05 


2 


59 ± 0.15 





32 


2 


IS 


2.07 


KK16-01 


2158.83 


2265.25 


01:55:23.51 


27:57:29.12 


-5.27 ± 


05 


1 


05 ± 0.01 


2 


40 ± 


12 


3 


77 ± 0.21 





IS 


1 


20 


2.50 


KK27-03 


3744.75 


3274.22 


03:20:50.39 


-66:19:38.16 


-4.81 ± 


07 


1 


02 ± 0.01 


1 


96 ± 


05 


3 


08 ± 0.11 





00 


1 


57 


7.05 


N1311-01 


493.31 


415.26 


03:20:20.05 


-52:10:14.95 


-5.03 ± 


01 


1 


09 ± 0.01 


2 


06 ± 


06 


3 


24 ± 0.11 





OS 


3 


57 


2.95 


N1311-04 


325.69 , 


3554.20 


03:20:04.00 


-52:09:21.26 


-6.49 ± 


03 





82 ± 0.01 


3 


63 ± 


07 


3 


35 ± 0.11 





IS 


3 


05 


2.52 


N1311-05 


2219.00 


2235.00 


03:20:07.86 


-52:11:11.03 


-8.59 ± 


03 





91 ± 0.01 


2 


98 ± 


01 


2 


75 ± 0.06 





00 





27 


0.22 


N1311-06 


1997.27 


2526.72 


03:20:06.69 


-52:10:56.18 


-7.33 ± 


04 





80 ± 0.01 


8 


47 ± 


06 


7 


81 ± 0.11 





04 





47 


0.39 


N784-01 


358.61 , 


1688.27 


02:01:18.09 


28:48:42.92 


-7.53 ± 


05 





97 ± 0.01 


2 


54 ± 


00 


2 


34 ± 0.07 





17 


2 


35 


0.98 


N784-03 


1788.73 


2006.70 


02:01:16.42 


28:49:52.81 


-6.63 ± 


06 





88 ± 0.01 


3 


07 ± 


00 


2 


83 ± 0.07 





12 





59 


0.25 


N784-04 


3598.94 


3096.28 


02:01:16.57 


28:51:38.43 


-6.38 ± 


07 





73 ± 0.03 


10.59 ± 0.05 


9 


76 ± 0.12 





10 


2 


11 


0.88 


N784-09 


3802.58 


3811.71 


02:01:18.54 


28:52:05.14 


-6.61 ± 


06 





77 ± 0.02 


3 


75 ± 


03 


2 


81 ± 0.10 





09 


2 


82 


1.18 


U1281-02 


2147.49 


2346.95 


01:49:32.54 


32:35:25.26 


-7.58 ± 


04 





97 ± 0.01 


1 


82 ± 


02 


1 


36 ± 0.06 





11 





17 


0.12 


U3755-01 


1638.62 


2239.86 


07:13:51.95 


10:31:42.11 


-9.18 ± 


06 





80 ± 0.01 


1 


55 ± 


01 


1 


16 ± 0.18 





20 


1 


30 


0.96 


U3755-02 


2864.94 


2681.61 


07:13:50.61 


10:30:40.01 


-7.26 ± 


06 


1 


01 ± 0.01 


1 


73 ± 


02 


1 


30 ± 0.20 





12 


1 


16 


0.85 


U3755-03 


1497.29 


2691.75 


07:13:50.40 


10:31:48.33 


-6.58 ± 


06 





95 ± 0.01 


2 


53 ± 


07 


1 


90 ± 0.25 





19 


1 


72 


1.26 


U3755-04 


960.61 , 


2756.50 


07:13:50.11 


10:32:15.02 


— 8.73 ± 


06 


1 


02 ± 0.01 


1 


56 ± 


00 


1 


17 ± 0.18 





13 


2 


65 


1.95 


U3755-05 


2836.61 


2038.49 


07:13:52.78 


10:30:42.62 


— 6.64 ± 


08 


1 


06 ± 0.04 


4 


14 ± 


15 


3 


10 ± 0.33 





06 





06 


0.71 


U3755-06 


2003.00 


2200.00 


07:13:52.13 


10:31:23.97 


— 8.23 ± 


06 





89 ± 0.01 


3 


11 ± 


00 


2 


42 ± 0.18 





04 





66 


0.48 


U3755-07 


2236.00 


2229.53 


07:13:52.06 


10:31:12.28 


— 7.33 ± 


07 





86 ± 0.01 


3 


02 ± 


00 


2 


35 ± 0.18 





OS 





24 


0.18 


U3755-08 


2529.77 


2427.00 


07:13:51.43 


10:30:57.23 


— 8.25 ± 


07 





91 ± 0.01 


4 


53 ± 


00 


5 


98 ± 0.18 





If) 





40 


0.30 


U3974-01 


964.54 , 


1568.30 


07:41:59.49 


16:49:00.92 


— 7.20 ± 


01 





85 ± 0.01 


4 


93 ± 


05 


5 


41 ± 0.19 





16 


3 


06 


0.87 


U3974-02 


1910.55 


1800.60 


07:41:58.16 


16:48:16.12 


— 8.76 ± 


04 





97 ± 0.01 


1 


65 i 


01 


1 


50 ± 0.15 





14 


1 


57 


0.44 


U3974-03 


2419.29 


1739.74 


07:41:58.08 


16:47:50.52 


— 7.86 ± 


01 





94 ± 0.01 


2 


13 ± 


02 


1 


94 ± 0.16 





13 


1 


67 


0.47 


U3974-04 


2485.60 


2663.96 


07:41:54.87 


16:47:54.72 


-8.27 ± 


05 





99 ± 0.03 


1 


96 ± 


01 


1 


78 ± 0.15 





11 





63 


0.18 


U4115-01 


1217.34 


3760.26 


07:57:03.79 


14:22:41.01 


-7.61 ± 


05 





96 ± 0.01 


2 


16 ± 


02 


1 


96 ± 0.18 





06 


1 


39 


0.91 


U685-01 


2866.04 


2345.56 


01:07:26.18 


16:40:56.84 


-7.01 ± 


03 





94 ± 0.01 


3 


03 ± 


02 


5 


75 ± 0.06 





02 


1 


31 


1.88 


U685-03 


3015.40 


2596.32 


01:07:25.68 


16:40:44.19 


— 7.94 ± 


04 





98 ± 0.02 


4 


34 ± 


OS 


5 


70 ± 0.11 





07 


1 


22 


1.75 


U685-04 


2066.37 


2746.89 


01:07:23.60 


16:41:21.88 


— 8.63 ± 


03 





91 ± 0.01 


2 


13 ± 


00 


2 


80 ± 0.04 





13 





61 


0.88 


U685-05 


2993.85 


2926.33 


01:07:24.64 


16:40:37.05 


-7.82 ± 


03 





97 ± 0.01 


3 


75 ± 


02 


3 


59 ± 0.06 





04 


1 


01 


1.44 


TTPOC net 

U080-O0 


1991.53 


3155.90 


01:07:22.24 


16:41:15.14 


— 8.35 ± 


03 





nn i n ni 

90 ± 0.01 


2 


01 ± 


00 


1 


78 ± 0.04 





06 





25 


0.36 


E154-023-02 


500.587 


3732.66 


02:57:01.03 


-54:35:24.44 


-4.78 ± 


03 


1 


13 ± 0.01 


1 


87 ± 


02 


3.85 ± 0.03 





15 


2 


89 


1.09 


E154-023-03 


750.00 , 


3799.00 


02:57:01.06 


-54:35:11.54 


-4.39 ± 


03 


1 


28 ± 0.01 


2 


59 ± 


05 


5.32 ± 0.07 





00 


2 


71 


1.02 


KK17-01 


1821.97 


1410.22 


02:00:15.85 


28:50:42.59 


-4.58 ± 


08 


1 


22 ± 0.02 


2 


02 ± 


03 


3.56 ± 0.03 





03 


2 


12 


8.41 


KK27-02 


2665.38 


3563.55 


03:20:57.69 


-66:19:03.74 


-6.15 ± 


08 


1 


20 ± 0.02 


4 


55 ± 


01 


6.8 ± 0.01 





11 





58 


2.61 


N1311-02 


1570.35 


, 490.91 


03:20:17.96 


-52:11:05.37 


— 5.69 ± 


03 


1 


39 ± 0.01 


2 


13 ± 


02 


4.17 ± 0.02 





11 


2 


72 


2.24 


N784-07 


1398.57 


2354.43 


02:01:18.31 


28:49:44.61 


-5.71 ± 


06 


1 


22 ± 0.03 


3 


01 ± 


08 


5.62 ± 0.10 





00 





88 


0.37 


N784-08 


1020.28 


, 3049.7 


02:01:21.32 28:49:45.61 


-4.38 ± 


07 


1 


14 ± 0.04 


2 


79 ± 


08 


5.2 ± 0.10 





16 


1 


01 


0.67 


U3755-09 


2247.0 , 


2767.67 


07:13:50.24 10:31:10.73 


-6.62 ± 


06 


1 


22 ± 0.01 


5 


05 ± 


22 


13.44 ± 0.40 





12 





00 


0.66 


U4115-02 


2306.86 


, 890.00 


07:56:54.05 


14:23:40.30 


-5.72 ± 


05 


1 


20 ± 0.01 


3 


09 ± 


03 


8.55 ± 0.06 





16 


4 


71 


3.07 


TTyfl 1 t; n 1 } 
U4110-UO 


214 30 




07:56:59.20 14:21:53.09 


fi A A -V- n 
— 0.44 zn u 


uo 


^ 


97 4.(1/11 


3 


57 ± 


01 


9 


88 ± 0.01 


n 
u 


1 4 


3 


19 


9 98 


U41 15-04 


3215.00 


, 2926.0 


07:57:01.16 14:24:22.25 


-4.94 ± 


07 


1 


26 ± 0.04 


2 


17 ± 


01 


6.01 ± 0.02 





28 


2 


70 


1.76 


U4115-05 


925.00 , 


3822.00 


07:57:03.97 14:22:26.30 


-4.86 ± 


06 


1 


35 ± 0.02 


1 


54 ± 


05 


4.26 ± 0.10 





is 


1 


02 


1.25 


U8760-01 


2932.65 


3042.88 


13:50:50.73 38:01:48.27 


-4.80 ± 


07 


1 


28 ± 0.02 


3 


60 ± 


07 


4.16 ± 0.05 





22 





51 


0.89 


E115-021-01 


505.00 , 


1425.43 


02:37:54.24 


-61:18:45.39 


-4.78 ± 


03 


1 


48 ± 0.01 


2 


74 ± 


03 


4 


91 ± 0.04 





12 


2 


10 


1.52 


IC1959-01 


395.93 , 


1110.70 


03:33:20.36 


-50:23:35.17 


-5.43 ± 


05 


1 


46 ± 0.01 


3 


16 ± 


03 


6 


88 ± 0.04 





00 


3 


13 


2.61 


IC1959-03 


670.10 , 


1963.78 


03:33:15.70 


-50:23:39.77 


-5.42 ± 


05 


1 


54 ± 0.01 


2 


23 ± 


02 


4 


85 ± 0.04 





IS 


2 


20 


1.83 


KK27-01 


2377.27 


3224.14 


03:21:01.00 


-66:19:13.72 


-5.22 ± 


08 


1 


59 ± 0.01 


2 


84 ± 


(IS 


4 


25 ± 0.08 





04 





20 


0.88 


KK65-01 


2016.4 , 


2749.84 


07:42:39.75 


16:34:37.65 


-6.12 ± 


08 


1 


48 ± 0.02 


2 


75 ± 


01 


4 


32 ± 0.01 





21 





88 


0.87 


N1311-03 


2869.85 


2768.65 


03:20:04.05 


-52:11:34.39 


-6.47 ± 


03 


1 


66 ± 0.01 


2 


63 ± 


01 


5 


14 ± 0.01 





16 





86 


0.71 


N784-02 


3337.92 


1645.64 


02:01:12.29 


28:50:50.86 


-6.09 ± 


05 


1 


57 ± 0.01 


1 


83 ± 


01 


3 


41 ± 0.01 





13 


1 


80 


0.75 


U685-02 


3403.67 


2455.21 


01:07:26.77 16:40:30.71 


-4.66 ± 


03 


1 


74 ± 0.01 


1 


91 ± 


03 


3 


22 ± 0.04 





06 


1 


(is 


2.40 


U1281-01 


1439.38 


3836.35 


01:49:38.95 32:35:09.73 


-7.40 ± 


03 


1 


25 ± 0.01 


2 


75 ± 


02 


5 


05 ± 0.03 





10 


2 


13 


1.53 


U3974-05 


1729.83 


2874.48 


07:41:54.57 16:48:33.71 


-8.71 ± 


05 


1 


33 ± 0.02 


3 


08 ± 


01 


s 


88 ± 0.01 





11 


1 


07 


0.30 



TABLE 4 

Results of K-S tests for V — I, My, r h and <e between our bGCs and the old GCs in the LMC and our Galaxy. 





V — I 


bGCs 
My r h 




V -I 


LMC 

My r h 


1 


YH 


54% 


50% 1% 


1% 


25% 


34% 32% 


3% 


OH 


1% 


19% 86% 


3% 


1% 


78% 7% 


20% 


LMC 


73% 


66% 8% 


99% 









